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Abstract 

We have performed an exhaustive search for Gribov copies on the 
lattice and their possible dependence from finite temperature effects. 
We show that, for each value of lattice size, Gribov copies are dense 
in configuration space at low temperature but their density tend to 
lower when the temperature increases. We have investigated lattice 
sizes running from 16^ x 8 to 16^. 



Gribov [|^ has shown that Landau's gauge defined by the equation 

d-A = (1) 

does not fix the gauge uniquely, that is given a gauge orbit A, satisfying the 
Landau condition we can always find a gauge transformation g, such that 
also satisfies the equation (1). The gauge fixing definition can be sharpened 
further by following Zwanziger idea of requiring that, besides the Landau 
condition, all the eigenvalues of the Fadeev-Popov operator be positive. The 
domain defined by the condition d ■ A = 0, and |c? ■ D| > 0, is called the first 
Gribov horizon. 

Zwanziger was able to show that the Gribov horizon is a convex region 
and every gauge orbit crosses it at least once. From the point of view of 
weighting correctly configuration in the functional integral it would be im- 
portant to show that every orbit goes through the first Gribov horizon only 
once. This statement has to be given meaning within measure theory in 
function space, that is all orbit contributing to the functional integral with 
non zero measure must cross the first Gribov horizon only once. As of today 
this statement can't be proven in the continuum theory. By looking at the 
lattice formulation of the theory we can gain some insight. As many author 
have already shown, on the lattice we do have configuration with Gribov 
copies as well as configuration for which no Gribov's copy has been found. 
Since, every configuration partakes to the functional integral for every value 
of the coupling constant we have to discern whether they do with or without 
a measurable contribution. 

As we all know, the lattice theory presents two distinct phases in the 
physical temperature, and an interesting question is whether the measure 
on configurations admitting Gribov copies is temperature dependent and 
whether is affected by the temperature transition. Secondly and more inter- 
estingly, we would like to know the measure in orbit space of configurations 
admitting Gribov copies. 

Beside these theoretical aspects of the quantum field theory, there is a 
pragmatic reason behind this investigation. One of the most effective ways 
to increase the overlap between physical states and operators used to excite 
these states from the vacuum, is given by the smearing procedure [0]. This 
consist of constructing non local operators in a way that the physical extent 
of these operators is constant and does not go to zero with the lattice spacing, 



thus maintaining a good overlap between the operator and the desired phys- 
ical state. This procedure, though, is not gauge invariant, since any attempt 
to construct smeared gauge-invariant operators is very costly and not very 
effective. Therefore we fix the gauge, typically Landau or Coulomb, and we 
build our non local wave functions. This procedure would be questionable if 
we were to show that fixing to a given gauge does not uniquely determine 
the point in orbit space. 

We have studied the problem of Gribov copies in the four dimensional 
SU{3) theory with the standard Wilson action. We are in agreement with 
previous results (see 0] and 0) where it's found that the value of jS affects 
the density of Gribov copies. In effect we find that, while at low (3 these copies 
have an high density, at higher values they become much more rarefied, even 
if their density never goes to zero. The conjecture that Gribov copies would 
disappear at the critical temperature, does not seem to be true, and this 
qualitative behavior seems to be volume independent. 

In chapter one we briefly review Zwanziger's formulation of the Gribov 
problem and discuss the lattice implementation of the same, in chapter two 
we describe the numerical implementation of the gauge fixing algorithm and 
finally in chapter three we discuss our results. 



1 Landau's Gauge Condition on the Lattice 

Let us consider the functional 

F,[g] = J d'xTrAl{x)A^J{x), (2) 

where 

Af^ = g~'A^g + g'%g, (3) 

and 

g{x) = e--^^\ (4) 

is an SU (3) matrix and w is antihermitian and traceless. 
If we expand Fc in w{x) we have: 

Fc[g] = Fe[l] + J d^xTr[d ■ A{x)w{x)\ - J d^xTr [w\x)d^D^w{x)] + o{w^), 

(5) 



where 

D^w{x) = d^w{x) + [AfX^^,iu{x)], (6) 

is the covariant derivative and d-D^ is the Faddeev-Popov operator. Zwanziger's 
prescription to gauge fix consists in choosing as representative on the orbit 
the field A such that Fc attains a local minimum. It follows from () that 
Zwanziger's condition is equivalent io d ■ A — Q and d ■ D positive definite. 

This variational formulation of the gauge fixing condition can be easily 
generalized to the lattice. If we look at the functional: 



Fi[9\ = Re 



Y.TrUl{x) 



X,IJ. 



where 
and again 



m{x)^g{x)U^{x)g\x + f,), 



(7) 

(8) 
(9) 



g{x) = e-^(^). 

in the limit of a — > ( where a is the lattice spacing ) we have 

Fi[g]^F,[g] + o{A'). (10) 
Expanding in w{x) we can obtain the lattice definition oi d ■ A, that is 



d-Af^Tr 



T'^J2{u,{x)-Ul{x-f, 



(11) 



and the stationary condition will he d ■ Ai = 0. 

The lattice expression for the Faddccv-Popov operator is quite cumber- 
some, so we will not write it here but it is quite easy to obtain. 



2 The Gauge Fixing Algorithm and Search 
for Copies 

Our algorithm to gauge fix an SU (3) configuration is implemented on massive 
parallel machines, with a checkerboard updating of the gauge links. We start 



by observing that the functional Fi[g\ can be rewritten as a sum over even 
sites or odd sites only, for instance: 



Fi[g\^Re 



E Tr{Uf,{x)-rU%x-rtJi)) 



(12) 



and a similar expression for odd sites. 

Wc can now apply, alternatively, gauge transformations which are non 
trivial on even sites and the identity on odd sites and viceversa as we switch 
from the definition for Fi as a sum over even or odd sites. In this framework, 
each site contributes to the functional with an independent term and the 
functional dependence of each term from the gauge transformation is given 
by: 



Re 



g{x)Y,TT {u,{x) + Ul{x + f,)) 



(13) 



To implement our gauge fixing prescription we only need to produce an 
algorithm that at each step increases the value of the local sum alternatively 
on each checkerboard, since this implies that the functional Fi will increase 
monotonically. Since Fi is also bounded this will guarantee the convergence 
of the algorithm. There are many such algorithms and we have implemented 
a variety of them. 

We perform such a gauge ransformation g{x) first on a checkerboard and 
after on the other one. The iterations has been carried out until 



|9- < 10- 



-14 



(14) 



We, then, apply to the gauge fixed configuration U, a "large" random 
gauge transformation, then we repeat the gauge fixing procedure generating 
a new configuration U'. We are interested in whether U' is the same or not 
as U. In fact, since Fi does not distinguish between configuration that differ 
for a global gauge transformation, we are only interested whether U and U' 
differ for a global or a local gauge transformation. The simplest observable 
that can answer this question is 



^{U,U')^W[\Tr{u,ix)-%ix)) 



(15) 



Typically, when we do not find a copy, the value of A is less than 10 
otherwise, if we do find a copy it is of the order of 10~^ — 10~^. 



As an insurance against round-off errors we have kept track of the gauge 
transformation needed to fix the gauge and, at the end, apphed its inverse 
to the gauge fixed configuration. We have always been able to return to the 
original configuration within an accuracy of 10~^^ — 10~^^. 

3 Results 

We have generated configurations on lattices 16'^ x Nt, with A''^ = 8, 10, 14 and 
16, for a large set of /3 values. The configurations have been generated with 
the Local Hybrid Monte Carlo algorithm and were separated by 250 — 500 
trajectories. For all the /3 considered we checked, for a set of observables, 
that the autocorrelation has never been larger than 15. For each configu- 
ration we have fixed to the lattice Landau gauge, then, starting from the 
original configuration we have performed a random gauge transformation. 
After fixing this newly obtained configuration we have compared the two 
gauge fixed ones and determined whether we were in presence of a copy. 

The number of configuration tried and the number of copies is summarized 
in table 1. From the results reported on this table we could observe that, on 
all the lattices we have simulated, if at the lower values of /3 the percentage of 
Gribov copies over attempts is next to the 100%, at higher /3's it decreases to 
10 — 20%. This happens smoothly and we do not observe any sharp transition 
around the critical temperature. 
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Lattice 




No of Conf 


Found / Attempts 


16^ X 8 


5.8 


8 


40/40 




6.0 


8 


24/40 




6.2 


8 


6/40 




6.3 


8 


6/40 


16^ X 10 


6.0 


8 


30/40 




6.2 


8 


17/40 




6.3 


8 


6/40 


16^ X 14 


6.0 


8 


35/40 




6.2 


8 


10/40 




6.3 


8 


6/40 




6.4 


8 


3/40 


16=^ X 16 


6.0 


8 


35/40 




6.2 


8 


20/40 




6.3 


8 


19/40 




6.4 


8 


3/40 



Table 1: Summary of results 



